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Qh' Abstract 

'^U( ■ The gravitational impact inside and outside of a domain wall is studied in the context of f{R) 

^b ' gravity theory. The function f[R) is found which satisfies the stability conditions. Our results 

fvS , imply that the domain wall may cause topological inflation. 

^! 1 Introduction 



> 



When a symmetry is spontaneously broken by a phase transition, topological defects are produced 
(see e.g., [l] and [2] for a review). The generated defects are different depending upon the kind of 
broken symmetry. For example, in the case when a discrete symmetry Z2 breaks, a two dimensional 
planar defect is produced. Such defects are called domain walls (DW). If the spherical symmetry 0(3) 
is spontaneously broken to a rotational symmetry C/(l), then a point-like (zero-dimensional) defect 
^iZ . arises. It is a monopole. As the other defects, cosmic strings and textures are also known. They are 

|/~j ' the one- and three-dimensional topological defects, respectively. 

^O ■ It is assumed in the current cosmology that topological defects could be generated when the inflation 

^^ ' occurred. However, no topological defects have been discovered until now. If there are some defects 

f^ . somewhere in our universe, their density should be very small. Topological defects may also be the 

CO ' origin of the dark matter, the dark energy and the large-scale structure of our universe. Therefore, 

studies of topological defects are useful to understand the history of our universe. Because of that 

reasoning, topological defects were extensively studied in the past — see e.g., ref. [3] for the study 

t^ . of creation and evolution of topological defects, and refs. [3] and [S] about the gravitational fields 

%^ ' originated from the topological defects. In addition, an inflationary model by using topological defects 

C^ . (called the topological inflation) was proposed in [B] and 0. It is similar to the chaotic inflation [?], 

but the universe created by a topological defect continues to expand exponentially. Applications of 

DW to the dark energy were proposed in ^ . 

More recently, series of models describing inflation or dark energy were proposed by many authors. 
Usually, an explanation of accelerated expansion needs assumes a matter with a negative pressure. The 
slow-roll inflation ( iIU] and [IT]) assumes the existence of a new scalar fleld called the inflaton. There 
are also the theories describing the accelerated expansion of the universe by modifying the geometry of 
the universe. They are called the modified gravity theories. The simplest example is the f{K) gravity 
theory. 

The f{R) gravity is defined by replacing the R term in the Einstein-Hilbert action with f{R), 
where R is the Ricci scalar and f{R) is an arbitrary function in terms of R. At present, there are 
many inflationary and dark energy models by applying in f{R) gravity. The Starobinsky model |12| is 
well known as an inflationary model, and is still viable. A well known dark energy model is an ad-hoc 
model proposed by Appleby and Battye [13] [j. In addition to those models many other models were 



We studied applications of Starobinsky and Appleby-Battyc models in refs. 1141 and |15l respectively. 



proposed. There are many reviews by the different authors (see e.g., refs. 116] and |17l). 

Several studies applying the f{R) gravity to monopoles and DW were also suggested. In the context 
of General Relativity, monopoles and DW can be represented by a spherically and a plane symmetric 
metric. The de Sitter-Schwarzschild solution '18' is known as the exact static spherically symmetric 
solution in General Relativity. On the one hand, Multamaki [1^ found that the de Sitter-Schwarzschild 
solution is also valid in the f{R) gravity. Sebastiani [20] found an exact solution different from the de 
Sitter-Schwarzschild solution. On the other hand, as a plane symmetric solution in General Relativity, 
the Taub's solution [21 is to be noticed. In the context of the f{R) gravity, Sharif and Shamir [22| 
derived three constant-curvature solutions which correspond to the Taub's, anti-de Sitter and self- 
similar solutions, respectively. In addition, the gravitational fields of monopoles, cosmic strings and 
DW in the Starobinsky model were studied in ^23^. Most recently, Carames et al. .24i calculated the 
gravitational field of a monopole as the f{R) function. The gravitational field created by monopoles 
was suggested by Barriola and Vilenkin f35]. Carames et al. derived the corresponding f{R) function 
in the case considered in |25]. Next the natural question arises: how about DW? This paper is devoted 
to a calculation of the corresponding f{R) function in the case of DW. 

Most of the existing studies address DW under the assumption that the thickness of DW is negligible 
for simplicity. However, we consider here a thick DW and calculate the f{R) function outside and inside 
of DW. We employ the plane symmetric metric of Taub [H] . Widrow [55] studied the gravitational 
field of a thick DW in the context of General Relativity and calculated its gravitational field outside 
and inside the DW by using certain approximation. We essentially employ the same approximation 
[26] though in the context of f{R) gravity. 

This paper is organized as follows. To begin with, we give an outline about DW in General 
Relativity in Sec. [2] In Sec. [3] we briefiy review the metric formalism of the f{R) gravity theory. 
The equations of motion are derived in Sec. 21 Our main part is given by Sec. [5] Finally, Sec. [5] is 
conclusion. Throughout this paper we use the natural units, c — h — 1, and the space-time signature 
V^il' = (1,-1,-1,-1)- 

2 DW in General Relativity 

We recall how DW are described in General Relativity. As was mentioned in the Introduction, a DW 
is related to a discrete symmetry breaking. Such symmetry, for example, is present in the well-known 
Higgs-type potential. The corresponding Lagrangian is given by 

£ = i5^''9^09,(/)-l/(0), (2.1) 

where (/) is a scalar field and the potential V{(j)) is 

Vi^)^^icj^'~rjY. (2.2) 

In flat space-time y^^ = 77^,^, the classical DW solution for the Lagrangian ()2.ip is given by 

,^o = rytanh(^|-j, (2.3) 

where So 

measures the thickness of DW. The subscript "0" of the scalar field (po and the thickness 60 indicate 
those quantities "in flat space". As regards the Lagrangian above, its relevant equation of motion for 
the DW case takes the form 



(2.5) 



where the primes denote the derivatives with respect to z. Integrating this equation once yields 

i^o' = y{<l>^)- (2-6) 

The metric representing DW is plane symmetric j21 . Assuming that the DW lies on xy plane, the 
most general metric compatible with it is given by 

ds^ = A(t, z){dt^ - dz^) - B(i, z){dx^ + dy^). (2.7) 

In the case of a static DW, the A and B take the form [21] , 

/I = , - B = \ + Kz. (2.8) 

VI + Kz 

Actually, the space outside of a DW without thickness is flat except for z = [21]. Indeed, the Ricci 
tensor to be derived from eq. p.Sp vanishes. It is also worth noticing that, if DW has no thickness, 
it cannot be static, because such DW will collapse and finally become a singularity, say, a black hole. 
However, if DW has a thickness, a static DW is allowed [26 . 

We consider the case of a plane symmetric metric (|2.7I) . In this section we gave the equation of 
motion in flat space (12.51) . It changes slightly with the metric (|2.7p , as we are going to show in the next 
sections. We also compute all components of the energy-momentum tensor and of the Ricci tensor 
corresponding to the plane symmetric metric. However, next we briefly review the /(i?) gravity theory. 

3 j[R) gravity 

As was mentioned in the Introduction, the /(i?) gravity theory is the simplest modified gravity theory. 
It is the theory in which the i?-term in the Einstein-Hilbert action is replaced with an arbitrary function 
/(i?) in terms of the Ricci scalar R. The action is given by 



^Jd^xV^,f(R) + S^, (3.1) 



^^—2 

where k is the gravitational coupling constant, k^ = SttGn, Gn is the Newton's constant and S'm is the 
matter part. In the Einstein General Relativity one gets the Einstein equations as the gravitational 
equations of motion. In the f{R) gravity the equations of motion are more complicated, 

MR)R^. - \f{R)9^. + (.9,..n - V,,V,)/,,(i?) = K,%,,, (3.2) 

where fn denotes the derivative with respect to R, 






and T^i, is the energy-momentum tensor. 



r,. = ^|%. (3.4) 

Taking the trace of eq. p.2p we have 

fRiR)R - 2f{R) + SafRiR) = K^T. (3.5) 



In this expression General Relativity corresponds to f{R) — R and /ij(-R) = 1. Then eq. p.5p 
becomes R = —k^T, and R is directly related to the matter T. In the context of the f{R) gravity. 



fniR) 7^ const., so n//f(i?) 7^ 0. It means that the new propagating scalar degree of freedom to = fR{R) 
is present. The dynamics of this scalar field (often called "scalaron") is described by eq. p.Sp . 
From the expression of eq. (13.51) . we can find that f{R) is represented by fn{R), 

f{R) = I [fRiR)R + 3n/fl(i?) - n^T] . (3.6) 



Substituting it into eq. p.2p . we find 



2rTi _ dil-V 



fR{R)R^.. - V^V./fl(i?) - K^T^. = ^ [fR{R) - afR{R) - K^T] . (3.7) 

The right hand side of this equation does not depend upon any indices. So, defining the quantity 

with fixed indices, we can obtain the identity 

C^-C, = 0. (3.9) 

This identity corresponds to the gravitational field equations. We will use it in the following sections. 
The f{R) function is also subject to some stability constraints. They are called "the stability 
conditions," and are given by 

fniR) > 0, fRR{R) > 0, (3.10) 

in our notation. The former is the classical stability condition, which means that graviton is not a 
ghost. In addition, the latter is the quantum stability condition, which means that scalaron is not a 
tachyon. So, both together they guarantee that the system does not get an instability [57] ■ We take 
them into account when we calculate the f{R) function below. 

4 Equations of motion of DW 

We derive the equations of motion of DW and the gravitational field equations in the plane symmetric 
space-time in this section. The metric of this space-time is given by eq. (|2.7p . We suppose the DW 
to be static, in order to simplify our calculations, so that A = A{z) and B — B{z). In the context of 
General Relativity, a derivation of the equations with a plane symmetric metric was done by Widrow 
[26] . He studied the gravitational effects inside and outside of thick and non-static DW. The static 
plane symmetric solution in the f{R) gravity was considered by Sharif and Shamir [22j. 

We start with the equations of motion. It can be easily calculated from the Lagrangian (|2.ip . 

where we used 

g^'d^d.cf = -La^(y^g^''a,0). (4.2) 

Next, the energy-momentum tensor is defined by 

T^u = d^c^d,(j)~g^,C. (4.3) 

In fact, most of its components vanish. The non- vanishing components of the energy momentum tensor 
are 

T% = T\ - T^2 = ]^<i>'^ + T/(0), (4.4) 

TS = -i0'2 - V{<P). (4.5) 



The (0, 0)-component Tqq, gives the surface energy density of DW, 

} 

dzToo. (4.6) 



/o 
In flat space it becomes 



a^^-^VXv'. (4.7) 

Note that the constant K introduced in Sec. [5] can be related to a, 

K^'^. (4.8) 

Next we compute the Ricci tensor R^^^ of the metric (j2.7p . Similarly to the energy-momentum 
tensor, most of the components of R^\, vanish. The non- vanishing components of the Ricci tensor are 

n 1 /A" A'^ A'B'\ , , 

(4.10) 
(4.11) 

The R^i and R^2 are the same because of the plane symmetry on xy-plane. From eqs. (|4.9p - (|4.1ip 
we find that the Ricci scalar R is 

2I I ^4 I2- + -§- - -§: 



R\ 


= R\ 


1 B" 
2A B ' 






R\ 


1 
~ 2A ' 


M" A'^ 2B" 
\ A A^ ^ B 


Q,2 
52 


A'B' 
AB 



R^—A- — + ^r-^]- (4-12) 



Substituting the energy-momentum tensor and the Ricci tensor into the identity p.9p we obtain 
the gravitational field equations. In particular, Co — Ci = gives 

A" ^'2 B" A'B'\ , {A' B'\ ^, , , 

/« + -r - 17 /fi = 0' (4-13) 



^A A^ B AB } '"■ \A B 
and Co — C3 = gives, 



2A'B' 2B" B'^\ 2A' , „ ■, n . ^ 

15T /k + -r/fl - 2/fl - 2At'<^'' = 0. (4.14) 



AB B B^ J ■" A 

In addition we get from Ci — C3 = that 

A" A'2 B" S'2 A'B'\ „ /.4' S' 



/« + -r + U /« - 2/fl - 2'^'<^'' = 0. (4.15) 



A A"^ B 52 AB J'" \A B 

Although there are three gravitational field equations, in fact, one of them is derivable from another 
two. Hence, there are in fact only two independent equations, the third one is redundant. 

The full set of equations of motion and the gravitational equations cannot be solved exactly. So, 
we are join got use some assumptions and approximations in next section, in order to find a solution. 

5 DW in f{R) gravity 

In this section, we find an approximate solution to the differential equations obtained in the previous 
section. First, we suppose that the corrections of A and B is small, so that 

A{z)^l + a{z), B{z) = l + b{z), (5.1) 



where have assumed a and h of the order k? . Under this approximation we have 

A' R' 

This approximation is the same as that used in |24j . 
Next, we expand (j) along the hnes of [26] as 

= 00 + '«^cr , (5.3) 

where C(k^) is neghgible. Finally, as for the expression of fR{z), we assume 

fRiz) = l + ^oz, (5.4) 

where V'o is a positive "small" constant. There is no special meaning in the expression of /^(z). We 
have just taken it for simplicity. Given these approximations and assumptions, the gravitational field 
equations (j4.13p - (j4.15p yield the following two equations: 

a" = 0, (5.5) 

/3" + 2kVo' = 0, (5.6) 

where we have introduced 

a — a — b , p ~ a + . (5.7) 

We now solve these differential equation inside and outside of DW separately, in the case of a thick 
DW. Recall that 6o denotes the thickness of DW. We find that the inside and the outside of DW can 
be represented by z/Sq > 1 and z/5q < 1, respectively. We make stronger conditions, 

f » 1, f « 1, (5.8) 

Oo Oq 

keeping all our assumptions and approximations. 

Having obtained the correction terms a, b for the metric A, B from eqs. (|5.5p and (|5.6p . we 
calculate the corresponding f{R) function. First, let us summarize our procedure. We use the Ricci 
scalar R obtained in eq. (I4.12p . Applying the approximations of eq. (|5.ip to R we get 

i?«i(3/3"-a"). (5.9) 

We can now calculate R — R{z) outside and inside of DW by substituting a and /3, respectively. Next 
we can obtain z — z{R). By substituting z into eq. (j5.4p . we write down fn in terms of R. Finally, we 
integrate that equation and derive the f{R) functions. 

5.1 Outside DW 

First of all, we notice that eq. (15. 5p implies 

a — a — b = ciz + C2, (5.10) 

where ci, C2 are integration constants. They can be determined by comparing eq. ()5.10p with static 
DW metric A and B. A and B are given by eq. ()2.8p outside of DW. But we approximated A and 
B by eq. (15. ip . Those approximated A and B should be equal to the expression in eq. (|2.8p being 
expanded in terms of z. Therefore, a and b should be given by 

a^--Kz, b^Kz. (5.11) 



Now we find that 

cx = --K, C2=0, (5.12) 

or 

a = Kz. (5.13) 

Next, we consider the expression of j3. Since 4>q is given by eq. (|2.3p . substituting it into eq. (|5.6p 

yields 

2^2772 

The solution is 



P" ^ _fy_sech4(e-/^o). (5.14) 



P = -\kz+^kW 



^ - In (e^/*° + e-'/^°\ + (p//''« + e-^/^^V 



(5.15) 



We determined the values of the integration constants by evaluating A^ B in the small z limit. 

We are now in a position to derive the Ricci scalar R by substituting eq. (|5.5p and (|5.14l) into 
eq. dnH), 



R 



^^sech\e'/'"). (5.16) 



^0^ 



When z/Sf) 3> 1, we can approximate sech {e/^^°) sa 2'*e '^^/^o. Therefore, we obtain 

So. 2^ . 3kW ,.,^. 



Then eq. (jO)) yields 

4 ^^ <52(-i?) 
Taking into account the stability condition, Jr > 0, we also get 



fn{R) = l + ^in^-j,P^. (5.18) 



-i^< ' g'^ e4/^°^-° = i?„,ax. (5.19) 

"0 

If \R\ is larger than i?maxj the system becomes unstable. But the value e^^^''^^ is very large, since we 
have assumed earlier that V'o is very small. Therefore, the classical stability can be satisfied. We can 
also calculate fRR{R) by differentiating fuiR) with respect to i?, 

fnniR) = ^y (5.20) 

This is positive definite, and, hence, the quantum stability is fulfilled too. 
Thus the corresponding f{R) function is given by 

/(i?)=(^l-^ji?+^i?ln-^^p^+I.C. (5.21) 

We assume that ■0o is small. Then equation (|5.2ip shows the Einstein-Hilbert action plus a small 
correction. Though we have left the exponential term in the expression of /3", it is very small since 
z/5a ~> 1. Hence, if we take e'^^/*" sa 0, then j3" = 0. It leads to i? = 0. Therefore, our result is valid, 
since it is consistent with the known fact that the space-time outside of DW is fiat. The integration 
constant in eq. ()5.2ip may be interpreted as a cosmological constant. 



In the above derivation we merely considered the positions on a distance from DW. How about the 
vicinity of DW? To consider the environment near DW, we put 



z _, z 



where z is smah. Then eq. ()5.6p becomes 



/3" 



4e2 
<5o 



4622 



-50 



where e\ and 62 are constants, 



ei 



2e 



It yields 



Therefore, we have 



0.18, 62 



462/ 126ie2K"'?7^ 



0.76. 



/fl(i?)= 1 + ^0^0 + 



-00(^0 \ -0o<5o(-i?) 



462 / 126ie2K^?7 



2„2 " 



(5.22) 



(5.23) 



(5.24) 



(5.25) 



(5.26) 



Taking into account that V'o is small, /i? > is satisfied. Differentiating it with respect to i?, we find 
that 

is clearly positive. Thus, indeed, the system under consideration is stable. The /(-R) function can be 
derived as the Starobinsky-type expression (/(i?) ^ R + R"^) [T^] with 



/(i?) =(1 + ^0-50 + ^ I i?- 



V'O'^o 



.R^ + I.C. 



Ae2 J '" ' 2'ieie2K?rf 
The coefRcients here are 462 ~ 3.0 and 246i62 ~ 3.3. If they are approximated by 3, we get 



/(^)= l + -V^o'5o i? + 



3 '^^ "y ?iK?rf' 



R^ + l.C. 



(5.28) 



(5.29) 



5.2 Inside DW 

In this subsection, we focus on the inside of DW. In this case, cfi'^ can be expanded asu 



4.12 _ n 



51 



^2 I 1 51 



Equation (j5.6p can be rewritten as 

P" = 
Integrating it in terms of z we obtain /3 






-'o 



/3 = - 



2k V /^l 



1 z" 



--' -l2jj'+^-^" + ^^' 



(5.30) 



(5.31) 



(5.32) 



^ The simplest option is to leave only the constant term in the expansion. However, it leads to almost the same result 
as in the case of the outside DW. Hence, we take the expansion up to the second order. 



where C3 and C4 are the integration constants. The equation for a ()5.5|) has the same form as that of 
the outside of DW. Therefore, a can be written down in the form of eq. (15.101) , 

a = C5Z + ce, (5.33) 

with the two integration constants C5, cq. However, the new integration constants we have now cannot 
be determined since the metric inside of DW is unknown. But we can calculate the f{R) function 
corresponding inside of DW. For that purpose we need a" and /3", not a and /3. 
We can solve eq. (|5.31|) for z 

z±=±6oJl~-^{-R). (5.34) 

The stability conditions 



/i,(i?) = l±Vo<Jo\/l-T^(-^)>0, (5.35) 



■00(50 




-1/2 



/i?«(^)=±6^(l-7I^(-^)) >0' (5-36) 



tell us that z+ is valid but z_ is not, while the constraint 



3k 7] 



2^2 



-R< —^ (5.37) 

''o 

is required. The corresponding f{R) function is given by 

;,«,.«,?^(,__|_,_«,y'%,a (5,38) 

The constraint (|5.37p is satisfied when i? w 0. If DW exists in our universe, it is known that 
ry < 100 [keV] [9]. Hence, the result we obtained may apply to the DW at present. 

There are also the models which require a large rj value. One of such models is the model of 
topological inflation [6], [7], which is similar to chaotic inflation, but with the topological defects 
created in the early universe which are exponentially expanding. Those defects in the topological 
inflation may continue to expand somewhere in our universe. The required value of 77 should be 
77 ^ K^^ in this type of inflation [5J [71 [55]. For this rj, the stability conditions corresponding to our 
model (|5.38p can be satisfied. Then, if A ^ 1, the right hand side of eq. (|5.37|) is of the order k~^ and, 
therefore, eq. (|5.37p can be —R ^ Sn'^ri^/Sg. Hence, the second term of eq. (|5.38p can be expanded as 

f{R)= ^''''f'^\ il + i^oSo)R, (5.39) 

where the integration constant is absorbed into first term. 

5.3 Estimate for -ipo 

Of course, it is desirable to get the observational bound for our parameter ^o- According to CODATA 
[29] , the value of the Newton's constant Gn is given by 

Gn = 6.6738(80) x 10"" [mVkg • s^], (5.40) 

in the SI units. The solution outside DW (|5.2ip should duplicate the Einstein-Hilbert action in the 
low-curvature (small \R\) regime. When \R\ <C 1, the logarithm term can be dropped, so that 

HR) ^-(i-^)r. (5.41) 



Equation (|5.4ip means that the correction ■0o<^o is added to the Newton's constant. This correc- 
tion should be smaUer than 10"''^ in order to be undetectable. Therefore, we obtain the constraint, 
^PqSq < 10^"*. For instance, if A ^ 1, we have 

ipa<lO^^T]. (5.42) 

Hence, by using the present DW constraint 77 < 100 [keV], we derive V'o ^ 10 [eV] ^ 10^ [1/™]- In the 
case of the topological inflation 77 ^ k^^ ^ 10^^ [GeV], we get V'o ^ 10^^ [GeV] ^ 10^^ [1/™]- 

Of course, we get only the upper bounds for tpQ. If it takes a smaller value, the correction to the 
gravitational constant Gn will appear in the next digits. Therefore, we conclude that the value of "00 
needs to be measured more precisely. 

6 Conclusion 

In this paper we derived the gravitational field inside and outside of DW in terms of the f{R) functions. 
The gravitational field equations in the /(i?) gravity are more complicated than those in General 
Relativity, and they cannot be exactly solved. Our results are derived in the weak field approximation. 
Applying the stability conditions, we have chosen the adequate solutions. The solution outside DW is 
given by eq. (j5.2ip . and in the vicinity of DW is given by eq. (15.281) . Inside DW we got eq. (|5.38p . 

Our results satisfy the constraint to exist DW in the present universe. In addition, those can 
be applied to the models which require large ry value, like the topological inflation. In deriving our 
solutions, we have introduced the constant ^pQ which has the inverse length dimension and yields a 
correction to the Einstein gravity. We determined the upper bound on the absolute value of tpo from 
the known value of the gravitational constant Gn- The magnitude of this correction is thus very small. 
Therefore, more precise observations are necessary to detect it. Though ipQ is likely to correspond to 
the constant K present in the static plane symmetric metric, their precise relation still remains to be 
found. 

Our solution of inside DW should be useful for describing the topological inflation in detail. How- 
ever, we assumed the DW to be static in the present paper. In contrast, the DW in the topological 
inflation continues to expand. Therefore, we should take into account the time evolution of DW, which 
we hope to consider elsewhere. 
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